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X ■ Abstract. 

Analysis of the Stockholm group data on zero-pressure-gradient boundary flows, 
presented in the thesis of J. M. Osterlund (|http : //www . mesh . kth . se/ ^ j ens/ zpg/| ) 



is performed. The results of processing of all 70 mean velocity profiles are pre- 
sented. It is demonstrated that, properly processed, these data lead to a con- 
clusion opposite from that of the thesis and related papers: they confirm the 
Reynolds-number-dependent scaling law and disprove the conclusion that the 
flow in the intermediate ("overlap") region is Reynolds-number-independent. 
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1 Introduction 



Turbulence is the state of vortical fluid flow where velocity, pressure, and other flow fleld 
properties vary in time and space sharply and irregularly, and can be assumed to be random. 
The experimental investigation of individual realizations of such flows is impossible because 
the results are irreproducible: Experiments repeated under identical external conditions 
produce a different outcome. Therefore experimental investigations of turbulent flows can 
only provide their average properties. 

Less trivial and not always recognized is the following: What is of greatest interest in these 
experiments are intermediate asymptotic states of wider classes of flows, i.e., coherent, self- 
consistent fragments common to many different ffows. Other measurements reflect special 
properties of a set-up, which cannot be reproduced in other experiments. It may sometimes 
be useful to investigate a particular device (e.g. an atomizer) for an immediate practical 
purpose, but one should be cautious in transferring the results to different ffows. 

Typical examples of intermediate-asymptotic ffows are shear flows, where the ffow is ho- 
mogeneous in the direction of mean velocity which depends only on the coordinate perpen- 
dicular to the mean ffow. A well known example of shear ffows occurs in smooth cylindrical 
circular pipes far from the entrance and outlet. Another example is the zero-pressure gra- 
dient boundary layer above a smooth ffat plate far from its tip. In spite of their apparent 
simplicity experiments with such flows require high experimental culture and are expensive, 
and therefore relatively rare. When they are successful, like, for instance, the experiments of 
Nikuradze [^] with flows in smooth pipes in the range of Reynolds numbers up to 3.24 • 10®, 
performed 70 years ago under the guidance of L.Prandtl, they become milestones in turbu- 
lence studies. They are used to check theories based on special hypotheses valid for special 
classes of flows. Not enough is known now about the solutions of Navier-Stokes equations to 
avoid such hypotheses. 

The responsibility of experimentalists who perform such experiments and process the 
data is therefore very high. They are to be very careful in their conclusions. We showed, for 
instance p'"^] that the experiments of Princeton group (A.Smits, M.V.Zagarola) presented in 
the thesis of M.V.Zagarola ['^], which attempted to increase the range of Reynolds numbers 
achieved by Nikuradze by an order of magnitude, had a flaw. Starting from Re = 10®, well 
below the upper boundary of Nikuradze experiments, their data were influenced by roughness 
— insufficient polishing of the pipe walls. 

In the present work we analyze the experiments with zero-pressure-gradient boundary 
layers presented in the thesis of J.M.Osterlund p]. Like the earher theses of M.V.Zagarola 
[^] and M.H. Kites [®], the thesis of J.M.Osterlund presents the results of long-time work 
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of a group headed by a senior scientist (in this case A.V.Johansson), using a comphcated, 
expensive and unique facihty. The experimental data and their interpretation presented in 
this thesis might be accepted by some readers, as in the case of the thesis by Zagarola, 
without precautions. This was the motivation of our analysis of this thesis. 

The number of runs (70 measurements of mean velocity profiles) reported in the thesis p] 
is larger than in previously reported series, although the range of Reynolds numbers covered 
is not as wide; less than in the previous thesis of M.H. Kites f]: 2500 < Ree < 27,500. In 
the thesis [^] the authors make very definite statements: they claim that their experiments 
confirm the classical two-layer theory, in particular the Reynolds number-independent uni- 
versal logarithmic law, and, exhibiting no Reynolds number-dependence, disagree with the 
alternative theory based on the Reynolds number-dependent scaling law. 

J.M.Osterlund presented the data of 70 mean velocity measurements on the Internet 
(www . mesh . kth . se/~ j ens/ zpg/ ) . We present here the results of the processing of all these 
data. We demonstrate that, properly processed, these data lead to the opposite conclusion: 
they confirm the Reynolds number-dependent scaling law and disagree with the conclusion 
of Reynolds number-independence. 



2 Background 

According to the classical two-layer theory of wall-bounded turbulent shear flows at large 
Reynolds numbers, the distribution of average velocity u across the flow in the basic in- 
termediate region adjacent to the viscous sublayer is represented in the form of universal 
(Reynolds number-independent) von Karman-Prandtl logarithmic law 

(P=-\nT] + C. (2.1) 

Here we use classical Nikuradze-Schlichting et al. notations: 

= —, = (-)2^ r] = , (2.2) 

p V 

where r is the shear stress at the wall, y the distance from the wall; v and p are the fluid's 
kinematic viscosity and density. In the thesis more modern notations are used: tJ^ instead 
of 0, y'^ instead of 77. So, the von Karman-Prandtl universal law ( |2.1| ) is represented in the 
thesis in the form 

f7+ = - ln(y+) + 5 . (2.3) 
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Here k (von Karman constant) and B, according to the logic of the derivation, should be 
universal constants identical in all high quality experiments. It is known from the literature 
however that various experiments give substantially different values of these constants. Niku- 
radze [^] determined n = 0.417, B = 5.84, Monin and Yaglom [^] give the values k = 0.40, 
B = 5.1; Schlichting gives the values k = 0.40, B = 5.5. The difference is substantial, 
and for many years doubts have accumulated on the validity of the universal logarithmic 
law. 

In our papers (see e.g. [^],[^°]) the derivation of the universal logarithmic law was recon- 
sidered. It was shown that one of the basic assumptions is not quite correct, and on the 
basis of an alternative assumption, a different "scaling" (power) law was proposed: 

U+ = {Clin Re + C2)(i/+)"/'°^'^ (2.4) 

where the constants Ci, C2, c should be universal and Reynolds number- independent. Com- 
parison with the experimental data of Nikuradze has given the following values of the con- 
stants: 



so that the law ( p.4| ) is presented in the form 

+ 5a 



2„ "=21^ 

or, using the notation of Osterlund's thesis p] 

^+=(i=lni?e+0(y+)^. (2.7) 

Asymptotically, at i?e ^ 00, the specific choice of Re is of no importance: Re can be replaced 
by Re' = Const ■ Re, and the asymptotic form of (|2.7| ) will remain the same. However, for 
practical large, but not too large, values of Re, its actual expression is significant. It should 
be remembered that in our comparison with Nikuradze's data [^] we used his definition of the 
Reynolds number: Re = ud/v, where u is the mean flow velocity (bulk discharge rate divided 
by cross-section area), and d is the diameter of the pipe. Furthermore, it was recognized from 
the beginning [^^] that the law (|2.7|) is asymptotic (in the parameter l/\n.Re). It should be 
valid at large Re, but at lesser values higher terms in the expansion of the coefficients could 
be significant. 
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The experimental data of the Stockholm group (J.M.Osterlund, A.V.Johansson) are pre- 
sented in the form of graphs in the Iny^, plane suggested by the universal logarithmic 
law ( |2.3| ) (see Figure 5 on page 43 and Figures 13, 14 on pages 152-153 of the thesis [^]). 
They became available on the Internet (www . mesh . kth . se/ ~ j ens/ zpg/ ) in complete digital 
form. The data are presented in our table, parameterized by the authors by the parameter 
Reg = U9/v. Here U is the free stream velocity, 9 is the momentum thickness, a quantity 
measurable a posteriori, after the run is performed. 

3 Processing of the mean velocity data 

The first question to be answered was as follows: are the mean velocity data presented on 
the Internet compatible with some scaling law, not necessarily the law (|2.7| ). Therefore 
the data were plotted in the double-logarithmic coordinates {\gy^ ,\gU^). The result was 
instructive: the data outside the viscous sublayer form a characteristic shape of a broken 
line (see Figures 1-70). This shape is similar to the shape obtained for the experiments of 
the first group according to our classification [^^] where there was no influence neither of 
the external turbulence of the free stream nor of roughness. The Stockholm authors place 
the lower boundary of the intermediate region at y^ = 200. We found this value generally 
too high, and the standard value y^ = 70 ^ 100 seems to be more appropriate, however we 
marked the line y^ = 200 on all Figures 1-70. 

Thus, all experiments revealed two straight lines forming a broken line in the lgy~^, ^gU^ 
plane. These straight lines correspond to the scaling laws 

(I) Tj+ = A{y+r (3.1) 
(in the intermediate region adjacent to the viscous sublayer), and 

(II) W=Biy+r (3.2) 

(in the intermediate region adjacent to the free stream). The coefficients A,a,B,f] were 
obtained by standard statistical processing of data (see the table). The coefficients A, a of 
the straight line ( p.l| ) representing the scaling law for the mean velocity distribution in the 
basic intermediate region adjacent to the viscous sublayer are obviously Reynolds- number- 
dependent: For us this was not unexpected, because previous processing of all available 
experimental data for a much wider range of Reynolds number led us to the same conclusion 
(see paper [^^] which was known to the Stockholm group and referenced by them). Therefore 
we conclude that the validity of some Reynolds-number-dependent scaling law for mean 
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velocity distribution is unquestionably confirmed by the experiments of the Stockholm group 
as well. 

Note that because the Reynolds number range covered by the Stockholm group was 
not large, substantially less than the range covered by the other groups, in particular the 
Chicago group f], there would be a danger that they would not notice the Reynolds-number- 
dependence because the governing parameter is In Re, not Re itself. This is not the case: 
The i?e-dependence of the Stockholm experimental data is sufficiently strong to be revealed 
by proper processing. 

By the way, the authors could notice that their values k = 0.38 and (3 = 4.1 are substan- 
tially less than those presented in the literature as standard. However, if the law is universal 
Reynolds-number-independent, these parameters should be identical for all experiments of 
sufficient quality! 

The argument against the power law used by the authors (see paper [^^] reproduced in 
the thesis) is the following. They introduce the "diagnostic function" 

P _ 1/+ dm 
U+ dy+ 

Their statement, "The function F should be a constant in a region governed by a power law" 
is correct for a fixed Reynolds number. However, this is not true for the "diagnostic function 
averaged for KTH data" , which is shown in their Figure 6. 

We invite the reader to look at any of the Figures 1-70. It is clear that for each run 
F is a constant — look at the straight lines in the first intermediate region! However, this 
constant is different for different runs because the slopes of straight lines is Re-dependent! 
Indeed, the slope in the first region decays with growing Reynolds number. It is clear why 
F obtained by the authors is decreasing: the runs with larger Reynolds number and smaller 
slopes contribute more to larger . 

Now, when the validity of some Reynolds-number-dependent scaling law for the experi- 
ments of the Stockholm group is unquestionably established, we have to investigate whether 
this scaling law can be represented in the same form ( p.7| ) obtained for flows in pipes. But 
what is Re? We cannot assume it arbitrarily to be equal to Rcq. 

This effective Reynolds number Re should have the form Re = UA/u, where U is the 
free stream velocity, u the kinematic viscosity, and A a length scale which we cannot d priori 
identify with the momentum thickness 6, as there is no rationale for such identification. 
So, the basic question is whether one can find for each run such length scale A so that the 
scaling law ( |2.7] ) will be valid for the mean velocity distribution in the first intermediate 
region. A priori the very existence of such a length scale is under question, but if it does 
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exist, this means that the law ( pI7| ) is not specific to flows in pipes but can be a general law 
for wall-bounded shear flows at large Reynolds numbers. 

To answer this question one has to take the values A and a for each run, obtained by 
standard statistical processing of the experimental data in the flrst intermediate scaling 
region, and then calculate two values lni?ei and \nRe2 by solving two equations suggested 
by the scaling law {\i.7\), 

-^lnRei + - = A, — ^ =a. (3.3) 

v/3 2 ' 2lnRe2 ^ ^ 



If the values lni?ei and \nRe2 obtained by solving two different equations (|3.3| ) are close, 
i.e. if they coincide within experimental accuracy, it will mean that the unique length scale A 
can be determined so that the experimental scaling law in the region (|3.1|) , whose existence 
was proved before, coincides with the basic scaling law (|2.7|) . The table shows that indeed 
these values are close — for all Reg > 10, 000, the difference ln_Rei — \1aRe2 does not exceed 
3%. This allows one to introduce for large Reynolds numbers the effective Reynolds number 
Re according to the relation 



lni?ei = i(lni?ei + lni?e2), or Re = ^/Re^Re^ , (3.4) 

i.e., the geometric mean of Rei and Re2- This Reynolds number allows the deflnition of the 
effective length scale A, which plays a similar role in the scaling law for the boundary layer 
flow as does the pipe diameter for flow in pipes. We remind the reader that the momentum 
thickness is calculated by integration of the velocity proflle obtained experimentally: the 
calculation of the length scale on the basis of the measured velocity proflle is not more 
complicated. Naturally the ratio of two length scales 6/ A is different for different runs: both 
these quantities depend upon the details of flows, in particular they can depend in principle 
upon the distance between the tip of the plate and the point of observation. 

4 Checking universality 

The scaling law (|2.7| ) can be reduced to a universal form 

'4) = -\n[— =ln|/+ (4.1) 



a VV3 + 5a 

where a = jT^ie- This formula gives another way to check the applicability of the Reynolds- 
number-dependent scaling law (p.7|) in the intermediate region ( p.l|) . Indeed, according to 
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( [4.1| ), in the coordinates Iny"*", ip, all experimental points should collapse onto the bisectrix of 
the first quadrant. Figure 71 shows that all data for large Reynolds numbers {Ree > 15, 000, 
24 runs) presented on the Internet collapse onto the bisectrix with accuracy sufficient to give 
an additional confirmation to the Reynolds-number-dependent scaling law ( |2.7| ). For lesser 
values of Reg a systematic parallel shift is observed (Figures 72-74). Apparently in this case 
the choice of Re according to ( |3.4D is insufficient because at small Reynolds numbers the 
higher terms of the expansion could have some influence (see the paper by Radhakrishnan 
Srinivasan [^^]). 



5 Conclusion 

The thesis of J.Osterlund contains the following statements: 

1. (p. 22 of thesis), "The classical two layer theory was confirmed and constant values of 
the slope of the logarithmic overlap region (i.e. von Karman constant) and additive 
constants were found and estimated to k = 0.38, B = 4.1, and Bi = 3.6." 

2. (p. 29 of thesis), in fact the Introduction to the paper [14]: "Contrary to the conclusions 
of some earlier publications, careful analysis of the data reveals no significant Reynolds 
number dependence for the parameters describing the overlap region using the classical 
logarithmic relation." 

These statements are not correct. Our careful analysis of experimental data presented in 
the thesis performed in the present work leads to the opposite conclusions. 

1'. The results contradict the classical two-layer theory. The estimates of the 
constants obtained by the authors are substantially different from the standard 
values and this reason alone is enough to reject the assumption of universahty, 
the cornerstone of the classical theory. 

2'. In full agreement with our earlier publications, careful analysis of the data 
reveals significant Reynolds number dependence for the parameters describing 
the "overlap" region and confirms the Reynolds-number-dependent scaling law. 

The thesis of J.Osterlund was not the first investigation of this kind, many similar ex- 
perimental investigations were performed earlier covering a much larger range of Reynolds 
number. (In the thesis only one decade of Reg was covered: 2, 500 < Reg < 27, 000; in 
previous investigations, in particular in those refiected in the instructive review [^^], the 
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range 1, 000 < Ree < 200, 000 was covered). However, as we showed above, the accuracy of 
experimental data is sufficient to reveal Reynolds number dependence and correspondence 
of their data to the Reynolds-number-dependent scaling law ( |2.7| ) proposed by us. 

The experiments of Osterlund et al. allow an additional con&mation of the separation of 
the basic part of the flow into two self-similar regions (I) and (II) governed by the laws ( |3.1| ) 
and ( p.2| ). It is important that these two self-similar regions cover the whole boundary layer 
and not a small (1/6!) part of it where the universal law is expected by the authors to be 
valid. These experiments reveal a weak i?e-number dependence of the parameter (3 (Figure 
75): it decreased with growing Re. The data are not sufficient to come to a final decision, 
but they are in approximate agreement with the correlation 

/3 = ^ + 0.01. (5.1) 
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